Transport in simple networks described by integrable discrete nonlinear Schrodinger 

equation 
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We elucidate the case in which the Ablowitz-Ladik (AL) type discrete nonlinear Schrodinger equa- 
tion (NLSE) on simple networks (e.g., star graphs and tree graphs) becomes completely integrable 
just as in the case of a simple 1-dimensional (1-d) discrete chain. The strength of cubic nonlinearity 
is different from bond to bond, and networks are assumed to have at least two semi-infinite bonds 
with one of them working as an incoming bond. The present work is a nontrivial extension of our 
preceding one (Sobirov et al, Phys. Rev. E 81, 066602 (2010)) on the continuum NLSE to the 
discrete case. We find: (1) the solution on each bond is a part of the universal (bond-independent) 
AL soliton solution on the 1-d discrete chain, but is multiplied by the inverse of square root of 
bond- dependent nonlinearity; (2) nonlinearities at individual bonds around each vertex must satisfy 
a sum rule; (3) under findings (1) and (2), there exist an infinite number of constants of motion. As 
a practical issue, with use of AL soliton injected through the incoming bond, we obtain transmission 
probabilities inversely proportional to the strength of nonlinearity on the outgoing bonds. 

PACS numbers: 03.75.-b, 05.45.-a,05.60.Gg. 
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I. INTRODUCTION 

We shall investigate transport in networks with ver- 
tices and bonds which received a growing attention re- 
cently. The networks of practical importance are those 
of nonlinear waveguides and and optical fibers [IJ , double 
helix of DNA Q, Josephson junction arrays with Bose- 
Einstein Condensates ^ , vein networks in leaves 0, 0] , 
etc. 

Major theoretical concern so far, however, is limited to 
solving stationary states of the linear Schrodinger equa- 
tion, and to obtaining the energy spectra in closed net- 
works and transmission probabilities for open networks 
with semi- infinite leads |6l4ll | . Only a few studies treats 
the nonlinear Schrodinger equation on simple networks, 
which are still limited to the analysis of its stationary 
statejl2. 

With introduction of the nonlinearity to the time- 
dependent Schrodinger equation, the network becomes 
to provide a nice playground where one can see interest- 
ing soliton propagations and nonlinear dynamics through 
the network 
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■ 17[ , namely through an assembly of con- 
tinuum line segments connected at vertices. Although 
there exist important analytical studies on the semi- 
infinite and finite chains 
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works within a framework of nonlinear Schrodinger equa- 
tion(NLSE) [22, Q. The subject is difBcult due to the 
presence of vertices where the underlying chain should 
bifurcate or multi-furcate in general. 

Recently, with a suitable boundary condition at each 
vertex we developed an exact analytical treatment of soli- 
ton propagation through networks within a framework of 
NLSE [25]. Under an appropriate relationship among val- 
ues of nonlinearity at individual bonds, we found nonlin- 
ear dynamics of solitons with no reflection at the vertex. 
We also showed that an infinite number of constants of 
motion are available for NLSE on networks, namely the 
mapping of Zakharov-Shabat (ZS)'s scheme 24| to net- 
works was achieved. 

The extension of the scenario to the discrete NLSE 
(DNLSE) is far from being obvious. The standard 
DNLSE is not integrable and the integrable variant of 
the continuum nonlinear Schrodinger equation is the one 
proposed by Ablowitz and Ladik (AL)!!!, E^-iM]- AL 



we find little exact an- 
alytical treatment of soliton propagation through net- 



equation is the appropriate choice for the zeroth order 
approximation in studying the soliton dynamics pertur- 
batively in physically motivated models, such as an array 
of coupled optical waveguide s i2ftl a nd proton dynamics 
in hydrogen-bonded chains [3d . |3l|. The dynamics of 
intrinsic localized modes in nonlinear lattices can be ap- 
proximately described by AL equation 32| . Exciton sys- 
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terns with exchange and dipole-dipole interactions also 
reduce to AL equation in some hmiting cases 33| . The 
AL chain is integrable by means of the inverse scattering 
transform, and, together with the Toda lattice 3j|, con- 
stitutes a paradigm of the completely integrable lattice 
systems. 

AL equation for a field variable iJj on one-dimensional 
(1-d) chain is given by 



n+1 



0, 



(1) 



where 7 is the strength of nonlinear inter-site interac- 
tion and n denotes each lattice site on the chain. This 
equation can be obtained from the canonical equation 
of motion with use of the non-standard Poisson brack- 
ets. Equation ^ has an infinite number of inde penden t 
constants of motion and is completely integrable 
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However, there is an ambiguity in generalizing the AL 
model to networks: how can we define the inter-site inter- 
action at each vertex in order to see the infinite-number 
of constants of motion in networks? To keep the inte- 
grability of AL equation, should any rule hold for the 
strength of nonlinearity on bonds joining at each vertex? 
We shall resolve these questions in this paper and show 
how solitons of AL equation on networks will be mapped 
to that of AL equation on a 1-d chain. Once this mapping 
will be found, the integrability properties like the inverse 
scattering transform, Backlund transformation, etc, are 
automatically guaranteed, and will not be addressed in 
this paper. 

Below we shall show the completely integrable case of 
the AL equation on networks with strength of nonlinear- 
ity different from bond to bond. As a relevant issue, with 
use of reflectionless propagation of AL soliton through 
networks, we shall evaluate the transmission probabili- 
ties on the outgoing bonds. In Section II, using a pri- 
mary star graph (PSG) and defining a suitable equation 
of motion at the vertex, we shall address the norm and 
energy conservations. In Section III, we shall show a ba- 
sic idea of the soliton propagation along the branched 
chain, finding the connection formula at the vertex and 
the sum rule among the strengths of nonlinearity on the 
bonds, which guarantee the infinite number of constants 
of motions and complete integrability of the system un- 
der consideration. In Section IV, the cases of generalized 
star graphs and tree graphs are investigated. Section V is 
devoted to the investigation of an injection of AL soliton 
which bifurcates at the vertex and is decomposed into a 
pair of solitons with each propagating along the outgoing 



bonds, and we shall evaluate the transmission probabil- 
ities on the outgoing bonds. Summary and discussions 
are given in Section VI. 



II. NORM AND ENERGY CONSERVATIONS 
ON PRIMARY STAR GRAPH 

A. Ablowitz-Ladik(AL) equation on networks 

Let us consider an elementary branched chain (see 
FigUJi namely, a primary star graph (PSG) consisting 
of three semi-infinite bonds connected at the vertex O. 
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FIG. 1: Primary star graph. 3 semi-infinite chains Bi, B2 and 
B3 connected at a vertex O. 

We denote individual lattice sites as {k,n), where k = 
1, 2, and 3 mean the bond's number and n corresponds 
to a lattice site on each bond. For the first bond k = 
1, n is numbered as n G _Bi = {0, — 1, — 2, • • • }, where 
(1,0) means the branching point, i.e., the vertex. For 
the second (fc = 2) and third {k — 3) bonds, n varies as 
n G Bfc = {1, 2, 3, • • • }. (2, 1) and (3, 1) stands for the 
points nearest to the vertex. 

Discrete nonlinear Schroodinger equation (DNLSE) d 
la Ablowitz-Ladik (AL) is defined on each bond except 
for the vicinity of the vertex as 

#fe,« + {lpk,n+l + 1pk,n-l) (1 + lk\lpk,nf) = 0, (2) 

where {k,n) ^ {(1, 0), (2, 1), (3, 1)}. It should be noted 
that 7fe may be different among bonds. There is an am- 
biguity about the interaction around the vertex, which is 
resolved as follows: Let's first introduce Hamiltonian for 
PSG as 



H = - ^(l/'i_„Vl,"-Hl + C.C.) - ^^{i^*k^ni^k.n+l + C.C.), 
n=0 fe=2n=l 

(3) 

where at the virtual site (1, 1) we assume ^/'u = S2i/'2,i + 
53i/'3.i with appropriate coefficients S2 and S3. Then 
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Eq. ([2]) can be obtained by the equation of motion 

iipk,7i = {H,ipk.n} (4) 

at {k, n) ^ {(1, 0), (2, 1), (3, 1)}, with use of non-standard 
Poisson brackets 

{V-fc,™, V-fc',™} = {V-fc,™, i^l',n} = 0- (5) 

On the same footing as above, the equation of motions 
in Eq.Q at (1,0), (2,1) and (3,1) are given, respectively, 
as 

i^l,0 + (V'1,-1 + 52-02,1 + SsV'S,!) (l + 7l|^l,0p) = 0, 

(6) 

#/c,i + (s/cV'i.o + '0fc,2) (l + 7fc|V'fc,iP) =0, fc = 2,3. 

(7) 

The solution is assumed to satisfy the following condi- 
tions at infinity: 4'i,n — at n — —oo and 4'k,n — > at 
n — >■ -|-oo for /c = 2 and 3. 



B. Norm and energy conservations 

It is known that the norm conservation is one of 
the most important physical conditions in conservative 
systems. Since Eqs.([2]),([6|) and ([7]) are available from 
Hamilton's equation of motion with non-standard Pois- 
son brackets, the norm and energy conservations seem 
obvious. Below, however, we observe them explicitly. 
Extending the definition in the case of 1-d chain [23| . 
the norm for PSG is given as 



iV = 11^-11' 



3 

fc=l neBk 



— ^ In (l + 7fc|?/'fc,i 



Its time derivative is given by 
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k,n 



k=l neBk 



with 
Ak 



(8) 



(9) 



(10) 



1 + 7fc|V'fc,nP 

For (fc, n) ^ {(1, 0), (2, 1), (3, 1)} with use of Eq. © we 



have 



Ak,n = - {lpk^n1Pln+l " V'fc,„V'fc,n+l) 



where 

jk,n = I ('0fe,nV'fc,„+l - V'fe,„V'fe,n+l) (12) 

implies a local current. Firstly one observes 

^ ^ Ak,n = ii.o - 32^ - h,i, (13) 

k n 

where ^ means the summation over all sites on PSG 

k n 

except for the points (1, 0), (2, 1), (3, 1). 

Then, for (fc, n) = (1, 0), (2, 1), (3, 1), with use of Eqs. 
(jni and (O wc obtain 

^1,0 = sar- - V'^,oV'2a) 

+ S3- (V'l,oV'3,l - '/'l*,oV'3,l) - JhO (14) 



and 



Ak.i = jk,i - Sk- ('0i,o'0fc,i - ^i,o^k,i) (15) 



for fc = 2,3. Substituting Eqs.dTS]), dH]) and (US]) into 
Eq.®, we can see -^N = 0, i.e., the norm conservation. 
Therefore, for any choice of values S2 and S3 the norm 
conservation turns out to hold well. 

On the other hand, the energy for PSG is expressed in 
a symmetrical form as 



E = -2Re 



3 +00 



,n— — 1 k—2 n—1 

+ V'1*0(S2V'2,1+S3V'3,1)] • (16) 

To show that the energy is conservative, we see its time 
derivative 

— £; = -2Re^ + ^l.nV'l.n-n) 

n— — 1 
3 +00 

-2Re^ (^Vfe,„^fc,«+i + ■0fc,„Vfc,n-Hi) - 

fe=2 n=l 

-2Re |^Vl,0 (s2^2,l + 53^34) + V'l.O (S2V'2,1 + S3V'3,l) 

(17) 

With use of Eq. ^ we have 

—00 

n=-l 
—00 

=- Y (i+7iiv'i,nn 



Jk,n Jk,n—1 1 



(11) 



n=-l 

-V'u-iV'i.o, 



(18) 
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and 



n=l 
^ oo 



n=2 



(19) 



The first terms in the final expressions in Eqs. (|18p and 
(|19P are obviously pure-imaginary. Substituting Eqs. 
(0 and (dH) into Eq.([T71) and using Egs.® and 
we find: 



dt 



-E 



-2Re 



53^3 



Tpl.O (S2V'2,1 + S3V'3,l) + ^1,-1^1,0 + 
+ V'2,lV'2,2 + '/'3,lV'3,2 

2Re 



-(i+7ii^i,or) 



|S2V'2,1 + S3V'3,in 



fe=l 



(20) 



The last equality comes from the pure-imaginary nature 
of the expression in [•••]. Equation (PH)) is nothing but 
the energy conservation. 

Thus we have proved that the norm and energy are 
conserved for any choice of values S2 and S3. In general, 
however, other conservation rules do not hold. In the 
next sections we shall reveal a special case with appro- 
priate choice of S2 and S3 which guarantees an infinite 
number of conservation laws. 



III. COMPLETELY INTEGRABLE CASE 
A. Dynamics near branching point and sum rule 

Among many possible choices of S2 and S3 , there is one 
special case in which an infinite number of constants of 
motion can be found and DNLSE on PSG becomes com- 
pletely integrable. To investigate this case, we shall first 
add to each bond Bk {k = 1, 2, 3) a ghost-bond counter- 
part so that Bk + B'^ constitutes an ideal 1-d chain 
(see Fig. [2]). Then we suppose that the soliton solution 
of AL equation on PSG is given by 



1 



qk,n{t), fc = l,2,3 



(21) 
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t= 
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FIG. 2: Real bonds and real solitons (solid lines) and ghost 
bonds and ghost solitons (broken lines). 



where qk,n{t) are soliton solutions of DNLSE with unit 
nonlinearity on the ideal 1-d chain ( 23, 2^ 27[): 



iqn + iqn+l + 9n-l)(l + kn P) = 0, 



(22) 



with n being integers in (—00, +00). The solutions of 
Eq. (P^ may be different among three fictitious chains 
Sfc + B^ (fc = 1,2,3). 

Comparing Eqs. ([S]), © and (f^ . one can find at the 
vertex the following two equalities: 

^ <Zi,i(i) - ^g2,i(i) + ^93,iW, (23) 



^ qkfiit) = ^qi.oit), fc = 2,3. (24) 



Noting the spatio-temporal behavior of soliton solu- 
tions and to guarantee the equality in Eg. (1241) . qk.n{t) — 
Sk^J^qi.n{t) with fc = 2,3 should be satisfied for any 
time t and for any integer n, from which we obtain 

fik , ^ Fn 



Sk 



and 



1 or Sk 



qk,n{t) = qn{t), 



(fc = 2,3) (25) 



(26) 



namely, the solution qk,n{t) should be bond-independent. 
With use of Eqs. (US]) and ^ in Eq. ^ we have the 
sum rule among nonlinearity coefficients 71, 72 and 73: 



1 

71 



1 

72 



1 

73 



(27) 



Equations ([25|) . (|26l) and ([27| are the necessary and 
sufficient conditions to see Eqs.(|23|) and (f24|) . Thus, un- 
der the sum rule for nonlinearity coefficients in Ea.(|27p. 
the solution on PSG is given by a common (bond- 
independent) soliton solution of Eq. (f22|) multiplied by 
square root of the inverse nonlinearity coefficient. For 
example, the soliton incoming through the bond Bi is 



expected to smoothly bifurcate at the vertex and prop- 
agate through the bonds B2 and B^, as we shall see in 
Fig. m In the case that 71, 72 and 73 break the sum rule, 
we shall see a completely different nonlinear dynamics of 
solitons such as their reflection and emergence of radia- 
tion at the vertex, as will be shown in FigEl The initial 
value problem for such a case is outside the scope of the 
present work. 

We also note that the parameters S2 and S3 would cor- 
respond to ^ and ^ , respectively, in the preceding work 
[25| , although the derivations of the connection formula 
at the vertex are quite different between the continuum 
and discrete systems. In fact, S2 and S3 are introduced 
to define the inter-site interaction at the vertex and are 
not obtained from the norm and energy conservations, in 
contrast to the case of networks consisting of continuum 
segments 



B. An Infinite number of constants of motion 

It is well known that Ablowitz-Ladik (AL) equation on 
the 1-d chain has an infinite number of constants of mo- 
tion. Now we shall proceed to obtain an infinite number 
of constants of motion for general solutions of AL equa- 
tion on PSG. First of all, it should be noted that the 
solution on PSG can now be written as 

1 



{qn{t)\ne Bk}, fc = 1,2,3, (28) 



where q(t) stands for a general solution of AL equation 
((22)) and is restricted to each bonds Bk [k = 1, 2, 3). 

While we already proved the conservation of energy, 
we can generalize it to the general case: Without taking 
the complex conjugate, Eq. (jS]) can be explicitly written 
as 

—00 3 +00 

n— — 1 k—2 n— 1 

-V't,o('S2V'2,l+S3V'3,l)- (29) 

Substituting Ea. (l28)) into Eq. ([29l) . Z is rewritten as 

^ —00 3 +00 



71 



1 



k=2 n=l 



— %<ll 
71 



E 



Sk 



—L V7i7fc 



(30) 



k=2 



Using the value s^ in Eq. (P5| and the sum rule in Ea. (l?7l) . 
Eg. (15(11) reduces to the constant for the ideal 1-d chain 



30: 



+ 00 



(31) 



Therefore Z in Ea. (l29l) is a constant of motion, and its 
real and imaginary parts imply the energy and current, 
respectively. 

For other higher-order conservation rules, we can write 
them as 



3 +00 

+ E-E/-H{9«|nei?,}), (32) 

k=2 11=1 

with defined as expansion coefficients of the expres- 
sion (see Ablowitz & Ladik [2?! ) 

log(5i°^ + 9^^z' + + •••) = ft^z' + /i"^-' + • • • , 

(33) 

where (5!™'') are given by 

gi°)-l, gl^^ ^ Rn-lQn-2, 
(m) _ ^n-1 (m-1) _ V- (""0 (0 ^ - 9 3 4 



Rn-2 



1 = 1 



-qn+2- 



(34) 



(35) 



The relations (|34)) and ([35]) are obtained by solving Eq. 
(4.15) in 0, i.e., 

2 Rn+1 



n+2 



Rn 



'idn+l ~ 1) — 2 Rn+lQn, 



recursively with use of the expansion 



(36) 



(37) 



The right-hand side of Eq. (|32|) includes some unde- 
fined field variables in the ghost bond regions which must 
be defined as 



72 



'ip2 



V': 



3,n 



with 71 > 1, 



il^k.n ^ \l —i'k,n, fc = 2,3 with n < 0. (38) 
V 7fc 

The conservation laws in Eq. ([5^ follows from the na- 
ture of solutions (j28p and the sum rule for nonlinearity 
coefficients (l?7l) . 



For m = 1 we obtain current and energy conservation 
laws. At m > 2 we obtain higher order conservation laws. 
Some of higher-order constants of motion are as follows: 



71 



1 ^ 

' k=ln€Bk 

+ Y<nirk,n+l)') , (39) 

3 

C3=-Y. [(V'fc,„+2V'fc,n-l(l+7fc|V^fc,n+in + 



+ (■0fc,„+l)Vfc,«V'fc,n-l)(l +7fc|'0fc,nn 



(40) 



where field variables at lattice sites of the ghost bonds 
are defined in Eq. 




FIG. 3: Tree graph. Bi ~ (-oo, 0), Bn, B12 ~ (0,L), and 
Biij ~ (0, +00) with i, j = 1, 2, • • • . 



B\ = Biij...„i and number the lattice sites on this bonds 
as 1, 2, 3, • • • , N/^. On each branching point we assume 
the following conditions are hold 



IV. GENERALIZED STAR AND TREE GRAPHS 



7A 7Am 



(43) 



Now we proceed to explore soliton solutions of DNLSE 
on other types of graphs and explore the sum rule and 
conservation rules for solitons to propagate through these 
graphs. 

The above treatment on PSG is also true for more 
general star graphs consisting of N semi-infinite bonds 
connected at a single vertex. In such cases, the initial 
soliton at an incoming bond Bi splits into TV — 1 solitons 
in the remaining bonds, and the extended version of Eq. 
(EZl) is 



The solution is given by 



1 ^ 1 

- = E- 

71 7j 



J=2 

The solution is given by the equations 

1 



i'k,n{t) 



/Ik 



--<ln{t), 



(41) 



(42) 



where n = 0, —1, —2, • • • for the first bond [k ~ 1) and 
n — 1, 2, 3, • • • for other bonds (2 < fc < N). qn{t) is a 
soliton solution of Eq. ^22\i . Conservation laws for this 
graph can be obtained analogously as in the case of PSG. 

Another example of the graph for which the soliton 
solution of DNLSE can be obtained analytically is a 
tree graph in Fig. |3l Now we shall provide a soli- 
ton solution in this case. We denote bonds of graph as 



V'A,«(t) = —=qn+sAt), n e Ba- (44) 
VTfc 



Here sa number of lattice sites that soliton pass through 
from Bi to Ba. For the tree graph it is defined as 

si = sii = no, siij = no + Nu, 

SA = su-j.-im ^no + Nu^ h Nuj...i. (45) 



Below, applying the induction method we give a 
proof of conservation laws for soliton solutions of AL 
on tree graph. Let us denote the tree graph as 
G and assume the conservation laws to hold in G: 
n+SA{t)) — const. Then we con- 
struct an enlarged tree graph in the following way: First 
we choose the arbitrary point in the one of right-most 
semi-infinite chain B^ as a new branching point. Cut of 
semi-infinite part of this bond at the point iV$ and attach 
M semi-infinite bonds to this point. Namely the bond 
-B$ is now replaced by a finite bond i?$ connected with 
M semi-infinite bonds i?$m — {1,2,--- ,N^rn}, with 
TO = 1, 2, • • • , M. For the enlarged tree graph, constants 
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of motion are given by network PSG in Figll] we shall give conservative quan- 

tities for the solution in Eq. ()47|) under the sum rule in 
fn \qn+SAit))+'li^ Y /"''H^n+s^W) Eq.(l27|. First of aU, the uorm in Eq. dS]) tumed out to be 
SaGG-b^, neBA neBs reduced to the One for the 1-d chain with the nouliuearity 

, .r-^ ,, s constant 71 and thereby is given by 

{t)) 

m=l neBcirn 

^ = 2/3/71- (48) 

BAeG-Bis. ne Ba n=i Equation (|48| indicates that a narrow soliton has a larger 

M +00 norm than wider ones. As for the energy (E) and current 

n+s,i.+Af4, (0) {J): it is convenient to evaluate the combined quantity 

m=i n=i+iv^„. ^ Eq^ ([29]) ^ith use of sa and S3 given by Eq. ((25|). 

(-1 -1 \ rCi-i/ \\ In fact we have 

7* -2^7<E.,„I 2^ fi,'>{q„+s^+N^{t)) + const. 
m=l / T! = l+Ar^„ 



(46) 



E^-2Re{Z), J^2lm{Z). 



(49) 



It is clear that the final expression becomes constant 
under the sum rule (|43p. Thus, starting from PSG in 
Fig. [T] and repeating the above procedure, we can get 
the conservation rule for all tree graphs. 



Substituting Eq. (H7)) into Eq. and using the sum 
rule in Eq. (|27p . one obtains 



Z = —e"'" sinh/3 
71 



(50) 



V. TRANSMISSION PROBABILITIES 
AGAINST INJECTION OF A SINGLE SOLITON 

A relevant issue of the above discoveries is the trans- 
mission probability against injection of a single soliton. 
Here we calculate transmission probabilities for a single 
soliton which is incoming through a semi-infinite bond 
Bi and outgoing through the other semi-infinite bonds 
{Bi\l^l}. 

A single (bright) soliton on a graph, which takes the 
general form as in Eqs. (H^ and is described 

with use of AL soliton with 7=1 [26]: ipi.nit) lying on 
individual bonds Bi is given by 



and 



ipi,n{t) —Ji sinh/3sech[/3(n — 7io — wt)] 



neBi, / = 1,2,3,--- ,iV, 



(47) 



where w — — 2cosh/3 cos a, v = — (2//3)sinh/3sina, — tt < 
a<7r, 0</3<oo,0<(/)o <27r and uq are bond- 
independent parameters characterizing frequency, veloc- 
ity, wave number, inverse width of the soliton, initial 
phase and initial center of mass, respectively. Equation 
(|47l) indicates that a narrow soliton travels faster than 
wider ones with the same a. 

It should be noted that parameter values are common 
to each bond, except for {7;}. Choosing the simplest 



4 4 
E = cos a sinh/3, J = sin a sinh/3. (51) 

71 71 



As is seen from Eq. (|T7)) . the center of mass of the 
soliton (CMS) on each bond Bi is located a.t n — uq a.t t — 
0. However, lattice points on the individual semi-infinite 
bonds are defined on the limited interval. In particular, 
on outgoing bonds {Bi\l ^ 1}, their lattice points n are 
defined in the interval (1, -|-cx)). If no < 0, therefore, CMS 
on {Bi\l 7^ 1} is initially located outside of the real bonds. 
In such cases we call the soliton as a "ghost soliton". 
When CMS belongs to a real bond we use a term "real 
soliton". In Fig. [2] which corresponds to PSG in Figll] 
ghost solitons are plotted by broken curve while real ones 
by solid line. The soliton dynamics here is governed by a 
single characteristic time r = While for < t < r 

the soliton at Bi is a real one and those at B2 and B3 
are ghosts, for t < t the soliton at Bi is a ghost and 
those at B2 and B3 are real. At < = with —no ^ 1, the 
soliton lying on the bond Bi is exclusively responsible for 
the norm N. On the other hand, at t ^ 1, the solitons 
running through the bonds B2 and B3 are exclusively 
responsible for the norm. Therefore we can naturally 
define transmission probabilities at t — >■ +00. 

In general networks, transmission probability for an 
arbitrary sem-infinite bond Bi(l 7^ 1) at discrete time t 
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50 too 150 200 250 300 

TIME 

FIG. 4: Numerical result for time evolution of a soliton prop- 
agation through a vertex in PSG. Strength of nonlinearity at 
each bond are 71 = 1, 72 = 1-5, 73 = 3 satisfying the sum 
rule in Ea. (l27|l . Space distribution of wave function proba- 
bility is depicted in every time interval T = 10.0 with time 
used commonly in branches 2 and 3. Abscissa represents dis- 
crete lattice coordinates defined in Fig[T] Initial profile is 
Ablowitz-Ladik soliton in Eq. (|47P &t t = Q with parameters 
j3 — 0.1, a — 57r/4. Time difference in numerical iteration is 
At — 0.01 Bottom panel shows the time dependence of partial 
norms at each of 3 branches. 

that makes vt integers are defined as 



n—1 



+00 



In (1 + sinh^ /3sech^{P{n - no - vi))) 

^' n=l 



+00 



y — In ( 1 -f sinh^ Z^sech^ {iSn')) 
Nji .71 ' 



n' — l~nQ—vi 



(52) 



At vt — > +00, X]„'=i_„(,_t,t the last line in Eq. ([52)) 
tends to X]n'=-oo ^^"^ ^^^^ summation gives iV, i.e., the 
normalization of the soliton in the ideal 1-d chain with 



^ 0.8 



o 



o 

w 
w 

E 
w 

n3 




0.5 0.6 0.7 0.8 0.9 

7 J 72 



FIG. 5: Transmission probabilities (TP) as a function of ^ 
in PSG. Symbols and lines denote numerical and theoretical 
results, respectively. A solid line with • and a broken line 
with o correspond to T2 and T3, respectively. 



the nonlinearity coefficient 71 . Therefore 

71 



Ti = 



7i 



(53) 



Under the sum rules as in Eqs. (I?7l) . (pi and P5|) we 
have the unitarity condition 



N 



(54) 



/=2 



where the summation is taken over the semi-infinite 
bonds except for Bi. The result in Eq. ([53| means that 
the transmission probability is inversely proportional to 
the strength of nonlinearity in outgoing semi-infinite 
bonds. 

We have checked this result using a numerical sim- 
ulation of the discrete nonlinear Schrodinger equa- 
tion(DNLSE) on PSG in FigdJ We numerically iterated 
Eqs. ([2]), dH) and ^ with use of Eq.(l25|) and chose the 
initial profile in Ea.(|47p with 71 and uq = —150 as an in- 
coming soliton. Figure |4] shows the result in the case that 
the sum rule in Eq.([27|) is satisfied: the soliton starting 
at lattice point n = — 150 in the branch 1 enters the ver- 
tex at n = and is smoothly split into a pair of smaller 
solitons in the branches 2 and 3 with no reflection at the 
vertex. The velocity and width of the soliton have the 
definite value common to all bonds, and the squared peak 
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value of the soliton is proportional to 7fc , which are con- 
sistent with the result in Ea. (|T7l) . Bottom panel in Fig|3] 
shows the time dependence of partial norms at each of 
3 branches. With increasing time, the partial norms at 
branches 2 and 3 converge to the transmission probabil- 
ities in Eq. (|53| . 

In FiglS] transmission probabilities T2 and T3 are plot- 
ted as a function of ^ in the wider range of 71 and 71 
in the case satisfying the sum rule in Ea. (P7)) . We can 
confirm the linear law predicted in Eq.([53|). 



, 0.04 I 1 r 




-200 -150 -100 -50 

site number 
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site number 




. W)0(YYYYWYYVYYy\ 1 
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0.5 I 1 1 1 1 1 1 

0.45 - 

Total 

0.4 ^ 

^ 0.35 ■ I 

O 0-3 ■ \ 

^ 0.25 ■ 1 B2 

0.2 ■ 1/ ■ 

0.15 ■ \ B3 

0.05 ■ ,v B, 

' ' ■ ■ ' ' 

50 100 150 200 250 300 

TIME 

FIG. 6: Numerical result for time evolution of a soliton propa- 
gation through a vertex in the case of 71 — 0.5, 72 = 1.5, 73 — 
3, which breaks the sum rule. Initial profile and parame- 
ter values are the same as in Fig[4l In top panel, broken 
curves indicate a propagation of the reflected soliton. Bottom 
panel shows the time dependence of partial norms at each of 
3 branches. 

Figure [6] shows the result in the case that the sum rule 
is broken: ^ -I- ^ 7^ 1. In this case the soliton starting at 
lattice point n = — 150 in the branch I enters the vertex 
at n = 0, but is accompanied with both reflection and 
emergence of radiation at the vertex. It is very inter- 



esting that the velocity of the self-organized soliton have 
the definite value common to all bonds. In particular, 
the refiected soliton at the branch 1 has the same mag- 
nitude of velocity as that of the incident soliton. With 
increasing time, the partial norms at branches 1, 2 and 
3 would converge to the refiection (on Bi) and transmis- 
sion probabilities (on B2 and B^). For some other choice 
of 71,72 and 73 that breaks the sum rule (which is not 
shown here), the asymptotically(i ^ l)-equal velocity of 
solitons running on all three semi-infinite bonds can also 
be observed and provides an open question to be resolved 
in due course. 

VI. SUMMARY AND DISCUSSIONS 

We have derived conditions under which Ablowitz- 
Ladik (AL) type discrete nonlinear Schrodinger equation 
(DNLSE) on simple networks is mapped to the original 
one on the ideal 1-d chain and becomes completely inte- 
grable. Here the strength of cubic nonlinearity is different 
from bond to bond, and networks are assumed to have 
at least two semi-infinite bonds with one of them used as 
an incoming bond. Our findings are: (1) the solution on 
each bond is a part of the universal (bond-independent) 
soliton solution of the completely-integrable DNLSE on 
the 1-d chain, but is multiplied by the inverse of square 
root of bond-dependent nonlinearity; (2) the inverse non- 
linearity at an incoming bond should be equal to the sum 
of inverse nonlinearities at the remaining outgoing bonds; 
(3) with use of the above two findings, there exist an in- 
finite number of constants of motion. The parameters 
S2 and S3, which played an essential role in deriving the 
connection formula, are introduced to define the inter-site 
interaction at the vertex and are not obtained from the 
norm and energy conservations, in marked contrast to the 
case of networks consisting of continuum segments [25}. 
The argument on a branched chain or a primary star 
graph (PSG) is generalized to general star graphs and 
tree graphs by using the induction method. As a prac- 
tical issue, with use of AL soliton injected through the 
incoming bond, we obtain transmission probabilities in- 
versely proportional to the strength of nonlinearity on 
the outgoing bonds. 

Acknowledgments. K.N. is grateful to F. Abdullaev for 
suggesting a significance of extending our preceding work 
on the continuum NLSE to the discrete case. 
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